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By determining the quark momentum fractions of the octet baryons from Nf = 2+1 lattice 
simulations, we are able to predict the degree of charge symmetry violation in the parton distribution 
functions of the nucleon. This is of importance, not only as a probe of our understanding of the 
non-perturbative structure of the proton but also because such a violation constrains the accuracy 
of global fits to parton distribution functions and hence the accuracy with which, for example, cross 
sections at the LHC can be predicted. A violation of charge symmetry may also be critical in cases 
where symmetries are used to guide the search for physics beyond the Standard Model. 
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Charge symmetry is related to the invariance of the 
QCD Hamiltonian under rotations about the 2-axis in 
isospace, turning u quarks to d and protons to neutrons. 
Extensive studies in nuclear systems have shown that it 
is an excellent symmetry [1] , typically accurate to a frac- 
tion of a percent (e.g. m„ — m p ~ 0.1%). At the quark 
level it is, of course, very badly broken but this is hid- 
den by dynamical chiral symmetry breaking. There has 
been extensive theoretical work on the effect of the u — d 
mass difference on parton distribution functions (PDFs) , 
where charge symmetry (CS) implies [2/3 : 

u?(x, Q 2 ) = d n (x, Q 2 ), dP{x, Q 2 ) = u n {x, Q 2 ) . (1) 

Within the MIT bag model, Sather [4] and Rodionov et 
al. [5] found that charge symmetry violation (CSV) in the 
singly represented valence sector, Sd(x) = d p (x) —u n (x), 
could be as large as 5% in the intermediate to large range 
of Bjorken x. Furthermore, these authors also found that 
5u{x) = u p (x) — d n (x) was similar in magnitude but of 
opposite sign. 

Only recently has a global analysis of PDFs allowed for 
CSV, with Martin et al. [5] finding a best fit that is re- 
markably close to the predictions of Ref. |S] for both the 
magnitude and shape of 8d(x) and Su(x). Unfortunately, 
the errors on their result are currently too large to be of 
phenomenological use but at the larger end CSV could 
lead to considerable uncertainties in the predictions for 
some processes of interest at the LHC. The need for ur- 
gency in obtaining better constraints on CSV in PDFs 
has recently become apparent in connection with the 
search for physics beyond the Standard Model using neu- 
trino deep-inelastic scattering. Indeed, the level of CSV 



predicted in Refs. [3J [5] would reduce the 3<r discrepancy 
with the Standard Model reported by the NuTeV col- 
laboration [7] by at least one standard deviation [5] . 
It was argued by Londergan and Thomas that for the 
second moments, which are relevant to the NuTeV mea- 
surement, namely (xSd~~(x)} and (x5u~(x)} (where the 
superscript minus indicates a C-odd or valence distribu- 
tion function), the results had very little model depen- 
dence |10j . Further, future planned new-physics searches 
will benefit from improved constraints on CSV, such as 
the parity-violating deep inelastic scattering program at 
Jefferson Lab [TT] . 

In this Letter we report the first lattice QCD deter- 
mination of the CSV arising from the u — d mass differ- 
ence. Our results are deduced by studying the second 
moments of the parton distribution functions as we vary 
the light (degenerate u, d) and strange quark masses in a 
Nf = 2 + 1 lattice simulation. The sign and magnitude 
of the effect which we find are consistent both with the 
estimates based on the MIT bag model [10] and with the 
best fit global determination of Ref. [6]. However, the 
uncertainties in this work are considerably smaller than 
those derived from the global analysis. 

Because of valence quark normalisation, the first mo- 
ments of 5u~ (x) and 6d~ (x) must vanish. Hence the sec- 
ond moment (which we label Sq~ ) is the first place where 
CSV can be visible in the valence quark distributions, 

Su- = [ dxx{vP-{x) - d n -{x)) = (x) p u _ - (i)3_ ,(2) 
Jo 

Sd- = [ dxx{dP-{x) - u n -(x)) = (x) p d _ - {x) n u - .(3) 
Jo 
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m K [MeV] 


(x)*/(x)Z 


(x)f/{xf d 


{x)f/{x)l 




0.12083 
0.12090 
0.12095 
0.12100 
0.12104 


0.12104 
0.12090 
0.12080 
0.12070 
0.12062 


460(17) 
423(15) 
395(14) 
360(13) 
334(12) 


401(15) 
423(15) 
438(16) 
451(16) 
463(17) 


1.0263(51) 

1.0 
0.9888(44) 
0.9670(83) 
0.9631(94) 


0.960(12) 

1.0 
1.0344(70) 
1.059(14) 
1.082(18) 


0.993(23) 

1.0 
1.010(25) 
1.019(26) 
1.037(29) 


1.044(28) 

1.0 
0.985(24) 
0.953(29) 
0.940(30) 



TABLE I. Pion and kaon masses on 24^ x 48 lattices with lattice spacing, a — 0.083(3)fm [12] . where the error on the lattice 
spacing has been included in the errors for and vtik- The last four columns contain our results for ratios of the hyperon 
quark momentum fractions. 



As detailed below, these CSV momentum fractions are 
related to the hyperon moments by 



5u 
8d~ 



(4) 
(5) 



in the limit where the strange and light quarks have al- 
most equal mass. 

In the numerical calculation of these moments, our 
gauge field configurations have been generated with Nf = 
2 + 1 flavours of dynamical fermions, using the Symanzik 
improved gluon action and nonperturbatively O(a) im- 
proved Wilson fermions |13j . The quark masses are 
chosen by first finding the SU(3)fl aV our-symmetric point 
where flavour singlet quantities take on their physical val- 
ues and then varying the individual quark masses while 
keeping the singlet quark mass m q — (m u +md+m s )/3 = 
(2m; + m s )/3 constant |12) . Simulations are performed 
on lattice volumes of 24 3 x 48 with lattice spacing, a — 
0.083(3)fm. A summary of our dynamical configurations 
is given in Table [TJ More details regarding the tuning of 
our simulation parameters can be found in Ref. |12j . 

On the lattice, we compute moments of the quark dis- 
tribution functions, q(x) 

{x n ~X = I dxx n -\q B {x) + {-l) n q B {x)) , (6) 



where x is the fraction of the momentum of baryon B 
carried by the quarks, by calculating the matrix elements 
of local twist-2 operators 



{B{p)\{0\ 



2(x n - 1 ) B 



-Tr] 
(7) 



where O^-^" 



f-^qy^D ■ ■■ D q . 
In this paper we consider only the quark-line connected 
contributions to the second (n = 2) moment, (x) q , which 
means we only include the part of q B coming from quark- 
line connected backward moving quarks, the so-called Z- 
graphs. While the contributions from disconnected inser- 
tions are expected to be small, in the following analysis 
we will focus on differences of baryons and so these con- 
tributions will cancel in the SU(3)fl avour limit and should 
be negligible for small expansions around this limit, as 
considered here. 



We use the standard local operator O 
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The matrix element in Eq. \7\ 



is obtained on the lattice by considering the ratio: 



R(t,T,p) 



C 3 pt(t, T,p) 

C 2p t(t,p) 



El 

V 



¥ 2 



(8) 



where Ci v t and C-$ v t are lattice two and three-point func- 
tions, respectively, with total momentum, p, (in our sim- 

-* \ (x) 

ulation we consider only p = 0). The operator O q is in- 
serted into the three-point function, C 3pt (t, r,p) at time, 
r, between the baryon source located at time, t — 0, and 
sink at time, t. 

The operators used for determining the quark momen- 
tum fractions need to be renormalised, preferably us- 
ing a nonperturbative method such as RI'-MOM [TH - 
I16j . Here, however, we will only present results for ratios 
of quark momentum fractions so that the renormalisa- 
tion constants cancel and hence our results are scale and 
scheme independent. In Fig. [l] we present results for the 
ratio of the w(s)-quark momentum fraction of the 
baryon to the momentum fraction of the u in the proton. 
They are also given in Table |lj as a function of m^, nor- 
malised with the centre-of-mass of the pseudoscalar me- 
son octet, = \J{2m 2 K +ml)/3 = 411 MeV. Here we 
see the strong effect of the decrease (increase) in the light 
(strange) quark momentum fractions as we approach the 
physical point. In particular, we see that the heavier 
strange quark in the S° carries a larger momentum frac- 
tion than the up quark in the proton. We also notice that 
the up quark in the S + has a smaller momentum fraction 
than the up quark in the proton. This is a purely envi- 
ronmental effect since the only difference between these 
two measurements is the mass of the spectator quark (s 
in S + , d in p). This implies that the momentum fraction 
of the strange quark in the S should be larger than that 
of the down quark in the proton, which is exactly what 
we see in Fig. [2j 

To infer the level of CSV relevant to the nucleon, 
we only need to consider a small expansion about the 
SU(3)flavour symmetric point, for which linear flavour ex- 
pansions prove to work extremely well [12] . For instance, 
we can write 



Su = ms I — 



d{x)l , d(x) 



dr 



drrid 



0(m 2 5 ), 



(9) 



where ms = (rrid — m u ) and we have already made use of 
charge symmetry by equating d(x}^/dm c i = d(x)P/dm u 
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FIG. 1. Ratio of doubly represented quark momentum frac- 
tions, {x)u/{x}Z and {x)f/(x)Z as a function of m^/X%, 
where we have determined X v from the masses in Tab. [I] 
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FIG. 2. Ratio of singly represented quark momentum frac- 
tions, (x)f I {x) p d and {x)^/{x) p d as a function of /X%, where 
we have determined X n from the masses in Tab. [i] 



and d(x)^/dm u — d(x) p u / dm d . A similar expression 
holds for Sd. 

Near the SU(3)fl avour symmetric point, we note that 
the up quark in the proton is equivalent to an up quark 
in a E + or a strange quark in a S°, which we describe 
collectively as the "doubly-represented" quark |17j . 

The local derivatives required for Su can be obtained 
by varying the masses of the up and down quarks in- 
dependently. Within the present calculation, we note 
that the difference (x)f — (x)? measures precisely the 
variation of the doubly-represented quark matrix element 
with respect to the doubly-represented quark mass (while 
holding the singly-represented quark mass fixed) . Similar 



variations allow us to evaluate the other required deriva- 
tives, where we write 



d(x)l _ (x)f -{x)l d(x)P _ {x)t~{x)l 



d(x) p d 
dm u 



mi 



mi 



dm d 

d(x) p d 
dm d 



m s 



m l 



mi 



-,(io) 



(11) 



With these expressions and Eq. ([9]), we obtain the rele- 
vant combinations for our determination of CSV in the 
nucleon 



Su = m$- 



mi 



Sd = m$- 



mi 



(12) 



By invoking the Gell-Mann-Oakes-Renner relation 
and normalising to the total nucleon isovector quark mo- 
mentum fraction, we write 



(13) 
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(14) 



Written in this way, the fractional CSV terms are just 
the slopes of the curves shown in Fig. [3] (evaluated at 
the symmetry point) multiplied by the ratio ms/m q . By 
fitting the slopes, we obtain 



Su 



I u—d 



Sd 



I u—d 



mg 
m q 
m s 
m„ 



(-0.221 ± 0.054) 



(0.195 ± 0.025) 



(15) 
(16) 



Chiral perturbation theory yields the quark mass ra- 
tio ms/m q = 0.066(7) [18j and the isovector momentum 
fraction is experimentally determined to be 



■t-d 



0.158 at 4GeV . Substituting these values into Eqs. (15) 



and ( 16 ) yields the first lattice QCD estimates of the CSV 



momentum fractions 

Su = -0.0023(6), Sd = 0.0020(3). 



(17) 



The first observation we make is that these results are 
roughly equal in magnitude and have opposite sign. 
These values are slightly larger than, but within er- 
rors in agreement with, the phenomenological predic- 
tions of E], where within the MIT bag model (at 



a scale Q ~ 4GeV ) they found Su 



-0.0014 and 



Sd = 0.0015. They are also consistent with the best- 
fit values of the phenomenological analysis of MRST [3] , 



Su = —Sd 



-0.002 



+0.009 
-0.006 



(90% CL). 



While our work provides the first nonperturbative 
QCD result to give a clear indication of the sign and mag- 
nitude of the CSV in these moments, we point out that it 
is based on lattice simulations using a single volume and 
lattice spacing. To achieve a precise quantitative deter- 
mination will require a detailed study of the finite- volume 
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FIG. 3. The difference between the doubly and singly repre- 
sented quarks in the E and S as a function of the strange/light 
quark mass difference. We deduce Su and Sd, respectively, 



from the slopes of these curves (c.f. Eqs. (15 1 and (161) 



and discretisation effects which we pfan to address by ex- 
tending these calculations to larger volumes and a second 
lattice spacing. Additionally, it is well known that lattice 
results for the second moment of the iso-vector nucleon 
PDFs, (x) u -d, do not agree well with experiment (see e.g. 
[19]). Based on chiral perturbation theory it is expected 
that finite size effects and chiral corrections are poten- 
tially large [2"0H2"3"] , but this has so far not been confirmed 
by lattice calculations. This discrepancy may also be due 
to a mismatch of lattice nucleon matrix elements and per- 
turbative Wilson coefficients. However, what concerns 
us here are ratios of moments of PDFs, in which such 
effects cancel out. For example, we find (x)f l /(x)^ l « 2.3 
in good agreement with {x) v _/{x) P d _ = 2.40(6) found 
in . We are also encouraged that lattice results for 
the ratio ( x )(Au-Ad) agree well with experiment 

[25] . Lastly, we have estimated the CSV associated only 
with the u — d mass difference. It is important to also 
find a method to investigate the CSV induced by elec- 
tromagnetic effects which is expected [S[ [55] to be of a 
similar size. The determination of this effect is, however, 
a separate calculation which will have no impact on our 
result. 

In summary, we have performed the first lattice deter- 
minations of the quark momentum fractions of the hy- 
perons, £ and S in Nf = 2 + 1 lattice QCD. By examin- 
ing the SU(3)fl aV our-breaking effects in these momentum 
fractions, we are able to extract the first QCD determi- 
nation of the size and sign of charge-symmetry violations 
in the parton distribution functions in the nucleon, 5u 
and 6d. Although our lattice calculations are restricted 
to the second (n — 2) moment of the C-even quark distri- 
butions, our results for Su = -0.0023(6), Sd = 0.0020(3) 



are in excellent agreement with earlier phenomenological 
calculations [8] . 
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